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Abstract 

We show that in the large Nc hmit four-quark twist-4 distributions in the pion can be 
expressed in terms of twist-2 pion distribution amphtude. This allows us to compute 
the isospin-2 structure function of the pion F2^'^{xb) in the large Nc limit. The method 
can be easily applied to other mesons as well. 

1. Matrix elements of four-quark operators appear e.g. in twist-4 contributions to the struc- 
ture functions of deep-inelastic scattering [H, 0, or in the description of non-leptonic weak 
decays . In this paper we consider the specific example of isospin-2 pion structure functions 

F2(xb) 



2Ff . 



2 — 2 "T 2 ^-^ 2 • (1) 

In this particular isospin combination of structure functions the contribution of twist-2 op- 
erators cancels out exactly and twist-4 contributions are dominant. It is remarkable that in 
the large Nc limit the twist-4 distribution function of quarks in the pion can be reduced to 
a convolution of twist-2 pion distribution amplitudes. 

2. One can read the contribution of four-quark twist-4 operators to F2 from |^ (correcting 
obvious misprints) 



Xb 
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dxidxdx2 Xb 



X [[^{x2,xi,x) + A(?/2,2/i,a;)]Tg(+)(a;i,x,X2) 
- [^{x2,yi,x) + A(2/2,a;i,x)]Tg(g")(xi,x,X2)} 
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Here yi 



X 



x\, 2/2 
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X2. We use the function A from ref. 



A(xi,X2,X3) = --Im 
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The functions Tjf^^(xi, x, X2) are defined as follows 



(3) 



r^^^^(xi,X,X2) = rqq(xi,X,X2) ± Tqg{xi,X,X2) 



(4) 



with Tgg, Tqq given by 



1 



Tgq{Xi, X, X2) 
Tqq{xi,X,X2) 



2P 



dzi dz 



2-K 2-K 2-K 



X exp \ ip 



XiZ^ + (x — Xi)z — (x — X2)Z2 I 

X (p|T{[^(o)n^^(z2-)] [^(z-)n^^(zr)]} b) . (5) 

Here g is the QCD coupling constant, = ^A^ are generators of the color group and 



7+(5e for Tjq 

Ibl^Qe for Tqq 



(6) 



where Qe is the quark electric charge matrix. The light-cone components of a vector are 
defined as = ^{V^ ± V'^). 

3. Generically the meson matrix elements of four-quark twist-4 operators are not related to 
matrix elements of two-quark operators. Here we show that in the limit of large number of 
colors Nc the matrix elements of four-quark operators can be expressed in terms of marix 
elements of two-quark operators. This observation allows us to express twist-4 contribution 
to ^2(0;) of the pion in terms of twist-2 distribution amplitudes. 

A simple analysis of leading in Nc Feynman diagrams P], corresponding to the meson 
matrix element 

(Mi|(^OiV')(^02V')|M2) (7) 

with matrices Oi being unity in the color subspace, shows that this matrix elements can be 
factorized in the limit of large number of colors as 



(Mi|(V^OiV^)(7/.02V^)|M2) (Mi|M2)(0|(V^OiV^)(^O2V^)|0) 

+ (Mi|(V>Oi7/>)|0) (0|(V^O2^)|M2) + (Mi|(V'O2^)|0) (0|(7/>OiV')|M2) . 
Let us define the connected part of the matrix element as follows 

connected 

= (Mi|(V^Oi^)(V^02V')|M2) - {M,\M2){0\{i>O,^){^O2n0) . 

If mesons Mi and M2 have different momenta or different internal quantum numbers then 
(M1IM2) = and there is no difference between the matrix elements and its connected 
part. But if |Mi) = IM2) then the nonconnected matrix element (Mi|(?/'Oi?/')('?/'02^/')|M2) 
is divergent and we must consider the connected part. Now the factorization takes the 
form 

ccmnected 



(8) 



(9) 



N, 



(Mil(V^OiV^)lO) (0|(V^O2V')|M2) + {M,\{i,02nO) mi'0MM2) . (10) 



We want to investigate the large behavior of the pion distribution functions. They 
are defined in terms of matrix elements 



(Mi|T{[^(4)rit^^(^i)] [V'(4)r2t''V^(^2)]} IM2) . 



(11) 
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We remind that 



A,B\ 



B 



Nr 



In the large A^^^ hmit this reduces to 



^^ii^jk ■ 



(12) 
(13) 

(14) 



Hence for any color singlet matrices Fi 2 we have 

(Mi|T{[^(zOrit^V^(^i)] yj{z'^)T2t''ij{z2)]}\M2) 



(15) 



We stress that here in ip{z'i) (S> ^"(^2) the tensor product refers to the spin flavor indices while 
the color indices are contracted. The trace Sp also refers only to spin-flavor indices. 

Applying the large Nc factorization formula ([ToD to the matrix element in the rhs of ( |15|) 
one obtains 



/ connected 



(Mi|T{[?/;(4)rit^^(zi)] [v5(4)r2t^v^(^2)]} IM2) 

-^Sp {{M,\T[ilj{z2)0ij{z[)] \0)Ti (0|t[^(zi)®V^(4)] |M2)r2} 

-^Sp{(0|T[V'(z2)®V'(4)] m)T,{M^\T[iP{z^)^^{z',)] 10)12} . (16) 

Now we apply the general large Nc eq. (|16D to matrix elements entering in definitions of 
twist-4 parton distributions (B|). We obtain the following result for these functions 



,X,X2} 



Tqg{xi,X, X2) + Tgq{xi,X, X2] 



i?2 

-fl'^-^ [^(^1 - 2^2 - l)07r(a;i - x)0^(Xi) + S{X2 -Xi- l)(l>n{x2 - X)4>n{x2)] 



4 
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(17) 



Here (pni^) is the pion twist- 2 distribution amplitude defined as 



1 



dX 



2n 



Here ~ 93 MeV is the pion decay constant and a, b are isospin indices. Also it is easy to 
see from definition eqs. (|,^ that at leading order of the 1/iVc expansion 



T^q \xi,x,X2) = 0. 



(19) 
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We see that in the large Nc hmit the meson twist-4 four-quark distributions are expressible 
in terms of twist-2 meson distribution amplitudes. Now we are in a position to compute the 
twist-4 part of F2{xb)- We insert expressions ( p!7||l9D into eq. @ 



-ituijst— 4 / 



[Xb, 



X 



-^2 P2 



207r(a;i?) / dxi4>^{xi) — 

J Xi 



- / dxidyi(t)^{yi)(t)^{xi)6{xi - yi - xb) 



Xi 

1 



+ 



1 



l-yi)yi (l-xi)xi 
From this equation we can easily compute pion isospin-2 structure function 

^0 



:^twist—A 



+ 



:^twist—A 



(20) 



(21) 



Here we have taken into account that the twist-2 contributions are cancelled in the rhs of 
eq. (|1|). The isospin factor of eq. (|1^) yelds 
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a3 rbS 
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a3 rb3 
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(22) 



Therefore the final expression for F2 (xb) has the form 



l — Xr 



dz(f)^{z) (f)^{xB + z) 
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z{l-z) 

1 



(23) 



_z{l - z) {xb + z){l - Xb - z) 
4. Using the large relation (^) one can see that our result obeys the normalization 

1 



constraint 



dxB ^{xb) 
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F2 
2_JL. 

2g2 



(24) 



independently of the shape of pion distribution amplitude (Ptt{x). We see that the first Mellin 
moment of F2^'^{xb) in the large Nc limit is given solely in terms of the pion decay constant 
and not sensitive to the details of internal structure of the pion. Recently this moment 
was computed in the lattice QCD M with the results 



1 

J dxB F^=\xb) = 0.27(10) g 



2FI 



(25) 



This calculation gives about two times smaller result for the first Mellin moment of F2 ^(a^s) 
than our model independent large Nc result. Taking into account various approximations 
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standing behind our result [leading order of the l/N^ expansion] and the lattice result 
( p5D [quenched approximation, extrapolation in light quark masses etc.] we find that the 
agreement between the two values is quite reasonable. Note that in contrast to the euclidean 
lattice calculations which can access the Mellin moments of parton distribution our approach 
allows us to compute the shape of the twist-4 distributions in mesons. 

Let us now compute the third Mellin moment of -F/"^(xb). This moment in the large Nc 
limit is sensitive to the form of the pion distribution amplitude 07r(x). To demonstrate this 
we expand the pion distribution amplitude in Gegenbauer series 



07r(a^) = 6x(l — x) 



i + E«nC'3/2(2x-i: 



n=2 



(26) 



where C^/^(2x — 1) are Gegenbauer polynomials. It is a matter of a simple calculation to 
see that the third Mellin moment of F/"^(xb) is sensitive only to the coefficient 02 in the 
Gegenbauer expansion of the pion distribution amplitude ( p6[ ) 



/ 



F^^ 3 



dxB xlFt\xB) = -{1 + as) . (27) 



Generically the A^th Mellin moment (for odd A^) can be expressed in terms of Gegenbauer 
coefficients up to order — 1. 

Using the general large N^. expression (^3[) we compute the twist-4 structure function 
F/"^(xs) for the asymptotic pion distribution amplitude 

0^(x) = 6x(l - x) . (28) 

The result (see Fig. |l|) is given by 

Ft\xB) = 6^ g'^xlil - xb){-1 + 2xb + 2x1) ■ (29) 

Q 

We see that function F^'^'^^xb) is concentrated at relatively large values of xb- Let us 
note that the isospin-2 pion structure function should have at least one zero due to the 
general sum rule 

f'^Ft\xB) = 0, (30) 
Jo Xb 

which directly follows from eq. (Q) and holds already at finite Nc- 

We have demonstrated that in the large Nc limit the higher twist quark-antiquark corre- 
lations in mesons are expressible completely in terms of their twist-2 distribution amplitudes. 
This model-independent result provides an illustration of general relations between higher 
twist structure functions and hadron wave functions as recently discussed by S.J. Brodsky [||. 
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Figure 1: The result for ■^^F2~'^{xb) in the large Nc limit for the case of the asymptotic 
pion distribution amplitude. 
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